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In the present paper, we investigate the properties of the sound generated by rubbing two objects. 
It is clear that the sound is generated because of the rubbing between the contacting rough surfaces 
of the objects. A model is presented to account for the role played by the surface roughness. The 
results indicate that tonal features of the sound can be generated due to the finiteness of the rubbing 
surfaces. In addition, the analysis shows that with increasing rubbing speed, more and more high 
frequency tones can be excited and the frequency band gets broader and broader, a feature which 
appears to agree with our intuition or experience. 
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Rubbing objects is an everyday experience. People 
who have gone through cold winters may have rubbed 
hands to get warm. Rubbing also occurs in Nature. The 
earthquake rupture is one of the most familiar but disas- 
trous ones. A common observation is that rubbing can 
generate sound or noise. The sound generated by rub- 
bing hands is certainly common to nearly every one. The 
sound generation by ice floe rubbing in the Arctic ocean 
contributes significantly to the ocean ambient noiseP,Q|, 
and may also help sea animals in finding appropriate 
holes or breaking segments in ice to breath. The sound 
generated by rubbing scoops against frying pans, how- 
ever, is likely something most people would like to avoid. 

In the present paper, we wish to consider the sound 
generated by rubbing two finite objects. A model is pro- 
posed for the sound generation due to the roughness of 
the contacting surfaces of the two objects. The sound 
field is calculated and is analyzed for its relation to the 
roughness and the rubbing speed. It is shown that sound 
field can be expressed in terms of a series of eigen-mode 
excitations, leading to tonal features. It is suggested that 
when the rubbing speed is low, only low frequency modes 
are possibly excited. With increasing speed, more and 
more high frequency modes can be excited, making the 
frequency spectrum broader. These features appear to 
be in accordance with our intuition or daily experience. 

Consider the problem of two object rubbing. The con- 
ceptual layout is shown in Fig. ^ One of the surfaces is 
moving with a constant velocity along the positive y axis, 
while the other is assumed to be at rest. For simplicity, 
we assume that the two objects are identical; for differ- 
ent sizes, the effective rubbing surface will equal that of 
the smaller one. The shape of the objects is a cubic with 
thickness L, width W and height H along the x-, y-, and 
2-th axes respectively. 

When rubbing occurs, a shear displacement field will 
be generated at the surfaces and is revealed as the shear 
waves. Such shear waves can leak out through defects 
or radiate out at the boundaries, and transform into the 
sound we hear or record by machines, e. g. microphones. 

The governing equation for the shear waves inside the 




FIG. 1: Conceptual layout of the system. 



objects is 
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where c s is the shear speed of the objects. 
The boundary conditions are 
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where fi is the shear modulus of the material. 

The stress can be separated into two parts: the mean 
and the variation, i. e. F(t,y,z) — (F(t,y,z)) + 
AF(y, z,t); here (•} denotes the ensemble average. We 
assume that (F(t, y, z)) is independent of y, z. It can be 
shown that the averaged stress can only excite the unin- 
teresting zero-mode, and will thus be ignored. 
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By Fourier transformation, 

u(x,y,z,t) = J duje~ tut u(x,y, z,u). (4) 

Eq. 0} becomes 

(V 2 +k 2 s )u( Xl y,z,u) = 0, (5) 
I 



with F(ui,y,z) = ^ J dte luJt F(t,y, z). It is easy to ver- 
ify that (A mn ), caused by the averaged stress, is propor- 
tional to 5 n ,o5m,o- Therefore only zero modes are possible 
for constant stresses. We will not discuss this situation. 

The purpose to calculate the intensity of the shear field 
which is related to the relevant sound field. That is, we 
need calculate the correlation function 

D(x,y,z,t;x',y',z',t') = (u(x,y,z,t)u*(x',y',z',t')). 

From the above derivation, it is clear that the key is to 
find the correlation function of the fluctuation part of 
the stress at the surface, (AF(y,z,t)AF(y',z',t')). Ob- 
viously the fluctuation is caused by the roughness of the 
contacting surfaces. If we assume that the roughness is 

I 



with 



The general solution to Eq. JSJ with Eqs. @ and JSJl is 
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homogeneous and completely random, i. e. spatially un- 
corrected at different point when the system is at rest. 
When the surface is moving against each other along the 
y-axis, the spatial separation along this direction is cor- 
related at a later time determined as A/ v. This consid- 
eration leads to 

{AF{y,z,t)AF{y',z',t')) = S8{z-z')5{y-y'-v{t-t')), 

(10) 

where S is a strength factor. 

Applying Eq. (|10|) . the intensity field can be calcu- 
lated. The procedure is to substitute Eqs. 10 and © 
into Eq. then calculate Eq. by taking into ac- 
count of Eq. iJTUjl. We finally get 



u(x,y,z,uj)= ^2 A 
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cos (A; (x - L)) sm(fc m y) cos(/c n z), 

n— 0,m— 
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where Then A mn is solved as 
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and m, n are positive integers. 

The coefficient A mn is determined from 



^ k mn A m ^ n sin(fc r „„L) sm(k m y) cos(k n z) = -F(v, \ 

(7) 
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A mn = r x — u x t ~i r-r, T~\~\ dz l dyF(uj,y, z) sin(k m y) cos(k n z), 



r 



< u(x,y, z,t)u*(x,y, z,t) > = S I du cos 2 (fc n z) y^ cos(k mn (x - L)) sm(k m y) cos(k m , n (x ~ L)) sxn(k m ,y) x 

n m,m' 

En , — ((-l) m e ^ w - l) ( (-l) m ' e l v w - l) (11) 
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This equation can be rewritten as 



where 



Therefore the frequency spectrum of the sound intensity 
field generated by rubbing is 

V ( w ) = -^y2 cos2 (knZ)\Q n (x,y,cj)\ 2 , (14) 

n 

with Q n being given by Eq. (|13() . 

A few general features can be observed from the spec- 
tral formula. First, the strength factor S controls the 
overall sound intensity level. It is expected to depend 
on a few parameters of the surfaces, such as the fric- 
tion coefficients and the mechanical properties of the 
surfaces. Second, due to the factor sin(fc m „L) in the 
denominator of Q n , the resonance feature is expected 
to appear when sm(k mn L) = 0, leading to the phe- 
nomenon of tonal sounds. This also shows that the thick- 
ness mainly defines the resonance feature. Next, the 
rubbing speed v enters Q n in the form of u>/v in the 
denominator. For a fixed frequency, the decreasing of 
moving speed will reduced the strength of Q n . There- 
fore with decreasing speed, high frequency components 
tend to decay accordingly. Furthermore, since k mn = 
y/{u}/c s ) 2 — (rmr/W) 2 — (nir/H) 2 , the cut-off frequen- 
cies are determined by (iv/c s ) 2 — {mir/W) 2 — (mr/H) 2 > 
0. These features comply with the intuition or experi- 
ence. It is apparent that the present model can be ex- 
tended to consider other rubbing surfaces, by adjusting 
the correlation function of the surface stress. All these 
features deduced from the theory tend to support the ex- 
perimental observation of the sound generated by ice-floe 
rubbing 2J. We note that the experimental results have 
been interpreted by an alternative model Q. 

To quantify our discussion, let us consider an example. 
Suppose that W = H. L = 0.1W, and S is a constant. 
Two moving speeds are considered: v/c s — 1/600 and 
1/300. As an example, we compute the sound field lo- 
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cated at (L, W/2, 0). The frequency spectra for the two 
moving speed are presented in Fig. El The units are arbi- 
trary. Here it is clearly shown that there is indeed a cut- 
off frequency at about 200 for the lower speed curve. Also 
for the lower speed case, the discussed resonance feature 
does appear. For the higher moving speed, the spectrum 
is broader. But the resonance tends to be weaker. And 
the cut-off frequency is larger (not shown). Compared 
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FIG. 2: Frequency spectra for two rubbing speeds. 



to the lower speed case, the frequency band is obviously 
broader. All these support the above qualitative analysis. 

In summary, a model is established to investigate the 
sound generated by rubbing two objects. In the model, 
we have assumed that the contacting surfaces are rough 
and the roughness is described as completely random. 
The dependence of the frequency spectrum on the moving 
speed and sample geometries is discussed. The results 
seem to be in line with intuition or experience. 
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